Seismic prospecting method and device using simultaneous emission of seismic signals based on pseudo-random sequences

ABSTRACT

Seismic prospecting method and device using simultaneous emission, by several vibrators, of seismic signals obtained by phase modulating a periodic signal whose amplitudes and derivatives with respect to the amplitude time are equated to zero at the beginning and at the end of each period, by pseudo-random sequences.  
     The seismic signals reflected by the subsoil discontinuities in response to the periodic signals emitted are picked up by receivers coupled with the formation and recorded in an acquisition and recording system. The seismic signals emitted are formed either from elementary sequences whose length is at least equal to the product of the number of seismic souces vibrating simultaneously by the listening time, or from this elementary sequence extended, before and after, with parts whose length is at least equal to the listening time, the reception and recording of the signals reflected by the subsoil discontinuities in response to the signals emitted, and the processing of the recorded signals. The respective contributions of the various seismic sources are isolated by correlating the signals received and recorded either by signals constructed from said elementary sequence extended, before and after, with parts whose length is at least equal to the listening time, obtained by circular permutation, or of said elementary sequence.

FIELD OF THE INVENTION

[0001] The present invention relates to a method and to a device intended for land seismic prospecting by simultaneous emission in the ground of seismic signals emitted by several vibrators or groups of vibrators, these signals being obtained by coding a signal by pseudo-random sequences, and notably a periodic signal phase modulated by such sequences.

BACKGROUND OF THE INVENTION

[0002] There are well-known land seismic prospecting methods comprising transmission in the ground, for several seconds, of a periodic signal whose frequency varies continuously within a frequency band, reception by pickups of the signals reflected by underground reflectors and recording of the signals received. As a result of the emission time, the signals picked up are combinations of signals reflected by reflectors arranged at very different depths. The image of the various reflectors in the subsoil can only be recovered by processing the signals picked up, including correlation thereof with the signals emitted. The processing result is identical to that obtained by convolving the autocorrelation function of the signal emitted by the reflection coefficients of the various reflectors. The seismic trace obtained is the image of the interfaces between the various geologic layers midway between the emission and reception points. Such a method is for example described in patent U.S. Pat. No. 2,688,124.

[0003] This method however has some drawbacks. The autocorrelation function that is obtained in this case exhibits secondary peaks on either side of the main peak, whose amplitude is quite significant. Furthermore, a time interval at least equal to the loop travel time of the waves emitted down to the deepest reflector in the zone explored, referred to as << listening interval>>, has to be provided between two successive emission sequences so that the strong signals picked up at the beginning of the corresponding recording sequence cannot conceal the weaker signals coming from more distant reflectors, picked up at the end of the previous recording sequence. The necessary emission interruptions during a relatively long listening time interval have the effect of limiting the energy transmitted.

[0004] There is also another well-known method wherein various vibrators emit simultaneously signals with an identical frequency sweep, the lag between their emissions being at least equal to the listening time.

[0005] Pseudo-random binary sequences

[0006] Other known land seismic prospecting methods use a technique that is well-known in the field of communications and radars. They comprise using periodic sources that transmit signals obtained by modulating the phase of a periodic carrier signal by a binary signal or pseudo-random code consisting of a sequence of elements that can take on two logical values 0 or 1. As shown in FIG. 1, the order of succession of these values is selected so as to have a random character.

[0007] When, for such a code, any sequence of n successive bits (n integer) can be repeated identically only after a sequence of (2^(n)−1) bits, the sequence is referred to as <<Binary Sequence of Maximum Length>> (BSML). These pseudo-random binary sequences can be used for phase modulation of a periodic signal, each element of the sequence being associated with a period of the signal, by keeping or by inverting the sign according to whether it is a <<1>> or a <<0>>. The term <<elementary sequence>> designates a portion of the periodic signal modulated by a binary sequence of (2^(n)−1) terms, of length (2^(n)−1).Tc, where Tc is the period of the carrier signal.

[0008] Correlation, by an elementary sequence, of a signal consisting of the repetition of this elementary sequence modulated by the elementary sequence (FIG. 2b) gives (FIG. 2c) peaks spaced out (in time) by the length of sequence T_(s) and a minimum (or even zero) level between the peaks or, more exactly, over the length of the sequence minus the period of the carrier T_(c). The ratio of the correlation peak to the correlation noise is equal to the number of terms of the sequence.

[0009] Such methods are for example described in patents U.S. Pat. Nos. 3,234,504; 3,264,606; 4,034,333 and 4,069,470.

[0010] U.S. Pat. No. 3,234,504 describes a seismic prospecting method wherein modulated signals without cusps, whose amplitude spectrum is centered on the zero frequency, are applied to a vibrator.

[0011] Patent FR-2,589,587 (U.S. Pat. No. 4,780,856) filed by the applicant describes a marine seismic prospecting method wherein the emission means comprise at least one vibrator towed by a boat progressing continuously along a seismic profile to be studied, emitting an uninterrupted series of sequences consisting each of a periodic carrier signal phase modulated by a pseudo-random binary coding signal of maximum length. The signals reflected by the discontinuities of the medium are correlated with the coded signals emitted so as to obtain correlation peaks at time intervals shorter than or at most equal to the period of repetition of the successive emission sequences. The periodic source can be the sole source, correlation being carried out between the signals received and alternately two sequences of signals emitted deducible from one another by a time lag shorter than the period of repetition of the sequences. A lag equal to the half-period of repetition of the sequences of signals emitted is for example selected.

[0012] It is also possible to use at least two periodic sources emitting simultaneously sequences of identical signals but with a time lag between them, and a correlation is established between the signals received, which correspond to the signals emitted simultaneously by the sources, and at least one sequence of coded signals, so as to alternately obtain correlation peaks corresponding to each periodic source.

[0013] The use, in land seismic prospecting, of vibrators emitting simultaneously or insufficiently separated in time has drawbacks linked with various factors: the autocorrelation noise, the harmonics and the slow waves.

[0014] Autocorrelation noise

[0015] The vibroseismic signal is compressed by correlating the signals recorded by the signal controlling the vibrator (or by a combination of the signals of the plate and mass accelerometers of the vibrator). The equivalent of the correlation of the series of reflection coefficients is thus obtained by autocorrelation of the signal emitted. The signal emitted is generally a frequency linear sweep whose amplitude spectrum has either a crenellated shape, or it is preferably a bell curve to reduce the amplitude of the bounces.

[0016] The question of autocorrelation noises arises for all the vibroseismic records. The bounces decrease with time as a function of 1/t substantially. For an isolated record, the autocorrelation bounces of the great values at the beginning of a trace are sufficiently attenuated when the weakest reflections at the end of a trace reappear. In the case of a frequency slip sweep, recording is semi-continuous and the bounces of the great values are found in front and behind, and they can interfere with the low values of the deep reflections of the previous shot if the slip time is insufficient.

[0017] Harmonics

[0018] For an isolated vibroseismic record, the harmonic distortion adds oscillations to the correlated signal. If the sweep is carried out from the low frequencies to the high frequencies, the oscillations due to the correlation of the harmonics by the control signal are precursors. Thus, except for the closest traces comprising the surface noise, the noises due to the correlation of the harmonics mix with an earlier, therefore in principle stronger signal. For continuous slip sweep type records, the noise due to the harmonics of the early arrivals of a shot can be superposed on the late and therefore weaker arrivals of the previous shot.

[0019] Slow waves

[0020] If the time interval between the start of two successive shots decreases, there is a risk the slowest waves of a shot, air waves and surface waves, may be found on the next shot. The sweeps being identical from one shot to the next, the air wave and the surface noises will be compressed similarly on the two shots.

[0021] Minimum correlation sequences

[0022] It is well-known to the man skilled in the art that there are, for each size of binary sequence of maximum length, minimum correlation pairs for which the ratio of the central peak to the greatest secondary peak is: $\frac{2^{n} - 1}{1 + 2^{{wholepart}{(\frac{n + 2}{2})}}}$

[0023] i.e. about 30 dB for n=11, corresponding to 2¹¹−1=2047 periods of the carrier signal (34 seconds for a 60 Hz carrier). There are connected sequences whose correlations are minimum two by two. The number of sequences of each set depends on the size of the sequences. For a sequence with 2047 elements, this number is 4.

[0024] Patent FR-2,818,753 filed by the applicant describes a method allowing to carry out prospecting or monitoring operations in an underground formation by elastic waves. It comprises emission, in the ground, of elementary sequences formed by coding of a signal by pseudo-random sequences, reception and recording of the signals reflected by the subsoil discontinuities in response to the signals emitted, and processing of the signals recorded by correlation with signals formed from the emitted signals.

[0025] The signals are emitted simultaneously by several seismic sources controlled each by a signal consisting of an elementary sequence whose length is at least equal to the product of the number of seismic sources vibrating simultaneously by the listening time, the respective contributions of the various seismic sources being separated by correlating the signals received and recorded by a cyclic signal portion constructed from said elementary sequence, extended, on the right and on the left (before and after), with parts whose length is at least equal to the listening time, obtained by circular permutation, this elementary sequence occupying the center of this cyclic signal portion.

[0026] The signals can also be emitted by several seismic sources controlled each by a signal consisting of an elementary sequence extended, before and after, with parts at least equal to the listening time, obtained by circular permutation, this elementary sequence occupying the center of this cyclic signal portion, the respective contributions of the various seismic sources being separated by correlating the signals received and recorded by signals constructed from said elementary sequence.

[0027] This method, which uses periodic signals obtained by phase modulation of a carrier signal by pseudo-random binary sequences, and this simultaneous recording of several seismic sources, allows to minimize the correlation noises and to prevent harmonics correlation noises. It is better suited to simultaneous recording than the frequency sweeps conventionally used in land seismic prospecting.

[0028] However, it has been observed that the phase modulation of a sinusoid by a pseudo-random code leads to cusps at code sign changes. The signal is more or less distorted upon emission. This distortion is expressed differently when passing between two like terms of the sequence and between two unlike terms. It follows therefrom that the various terms cannot be correctly summed, so that the correlation does not reach the desired performance.

SUMMARY OF THE INVENTION

[0029] The method according to the invention allows to carry out prospecting or monitoring operations in an underground formation by elastic waves. It comprises simultaneous emission, in the ground, of seismic waves by several seismic sources controlled each by a signal consisting either of an elementary sequence formed by coding of a signal by pseudo-random sequences, the length of this elementary sequence being at least equal to the product of the number of seismic sources vibrating simultaneously by the listening time, or respectively of this elementary sequence extended, before and after, with parts whose length is at least equal to the listening time, reception and recording of the signals reflected by the subsoil discontinuities in response to the signals emitted, and processing of the signals recorded, the respective contribution of each one of the various seismic sources being separated by correlating the signals received and recorded either by signals constructed from said elementary sequence, extended, before and after, with parts whose length is at least equal to the listening time, obtained by circular permutation, or respectively of said elementary sequence. It is characterized in that the signals applied to the various vibrators are formed from a periodic signal whose amplitude and derivatives with respect to the amplitude time are equated to zero at the beginning and at the end of each period, and such that its integral over a period is zero.

[0030] The integral over a period being zero, this prevents the frequency spectrum from being centered on the zero frequency (which of course cannot be emitted in practice) and therefore avoids the real signal deformations resulting therefrom.

[0031] According to an embodiment, a signal of frequency f, of the form sin 2π.ƒ.t.(1−cos 2π.ƒ.t), is used as the periodic signal.

[0032] With this additional condition applied to the periodic signal, the separation obtained of the contributions of the various seismic sources to the recorded signals is much better.

[0033] According to a preferred embodiment, the elementary sequences are formed by phase modulation of the periodic signal and these sequences are applied to vibrators.

[0034] The same control sequence is for example applied to all the vibrators, with a time lag and a circular permutation, the time lags between two sequences being at least equal to the listening time, and the length of the elementary control sequence being at least equal to the sum of the time lags and of the listening time.

[0035] According to another embodiment, the vibrators are controlled by a connected set of several minimum crosscorrelation sequences, each vibrator emitting a vibration with its own sequence independently of the other vibrators.

[0036] According to another embodiment, the previous two modes are combined by dividing the vibrators into several groups, and each group is assigned a pseudo-random sequence belonging to the same group of minimum crosscorrelation sequences, the various vibrators of the same group being controlled by the same elementary sequence with said time lags.

[0037] The device intended for prospecting or monitoring in a formation by elastic waves according to the invention comprises at least one group of m vibrators, a control unit suited to apply to the various vibrators elementary sequences formed by phase modulation of a periodic signal by pseudo-random sequences, seismic receivers coupled with the formation, a system intended for acquisition and recording of the seismic signals reflected by the subsoil discontinuities in response to the periodic signals emitted and a system for processing the seismic signals recorded, by correlation with part of the periodic signals emitted. The control unit comprises a periodic signal generator and a modulation set for forming elementary sequences phase modulated by a pseudo-random sequence, and connection means for applying simultaneously to the various vibrators of each group the signals produced by the modulation set.

[0038] According to a first embodiment, the modulation set comprises means for generating at least one pseudo-random control sequence, time lag means for forming at least one set of m elementary sequences with time lags between each one of them respectively, means intended for phase modulation of the periodic signals generated by said generator respectively by the m elementary sequences, generating m phase modulated periodic signals.

[0039] According to another embodiment, the modulation set comprises means for generating at least one pseudo-random control sequence, modulation means for modulating, by said sequence, the periodic signals generated by the generator, and time lag means for producing m modulated periodic signals with a time lag between each signal.

[0040] The device comprises for example p groups of vibrators, and the control unit is suited to generate p minimum crosscorrelation elementary sequences and to apply these sequences with a time lag to the vibrators of each group.

[0041] According to an embodiment of the device, the vibrators are installed permanently for seismic monitoring of an underground zone.

BRIEF DESCRIPTION OF THE FIGURES

[0042] Other features and advantages of the method and of the device according to the invention will be clear from reading the description hereafter, with reference to the accompanying drawings wherein:

[0043]FIGS. 1a to 1 b show examples of signals obtained by phase displacement,

[0044]FIG. 2 shows an example of a periodic signal whose time derivative is equated to zero at the beginning and at the end of each period,

[0045]FIGS. 3a to 3 c show examples of signals obtained by correlation by a sequence modified by a time lag,

[0046]FIGS. 4a to 4 e show various recorded signals corresponding to two phase modulated periodic sources that emit simultaneously,

[0047]FIG. 5 diagrammatically shows the device, and

[0048]FIG. 6 shows the control unit intended for the various vibrators.

DETAILED DESCRIPTION

[0049] We consider hereafter the formation of the control signals obtained by coding, by pseudo-random sequences, of a periodic signal a) that is equated to zero at the beginning and at the end of each period, b) whose derivative with respect to time is also equated to zero at the beginning and at the end of each period, and c) such that its integral over a period is zero.

[0050] Among all the possible signals having this property, a signal of the form 2π.ƒ.t.(1−cos 2π.ƒ.t) can be selected for example (FIG. 2).

[0051] Correlation by a sequence modified by a time lag

[0052] Let a sinusoidal signal of frequency ƒ_(c) (and of period $T_{c} = \frac{1}{f_{c}}$

[0053] phase modulated by

[0054] a BSML sequence (FIG. 1A). The length of the sequence S thus formed (FIG. 1B) is T_(s)=(2^(n)−1)T_(c), with n an integer. Sequence S is split up into two subsequences S₁, consisting of the 2^(n−1) first carrier periods of S, and S₂, consisting of the rest of the sequence, i.e. 2^(n−1)−1 carrier periods, S=(S₁ S₂) (FIG. 3a). Let S_(cyc) be the signal periodised from sequence S, i.e. consisting of cycles (S S S . . . ). Correlation of S_(cyc) by elementary sequence S gives (FIG. 3b) peaks spaced out by the length of the sequence T_(s), If we construct another BSML sequence by applying a time lag of 2^(n−1) terms and complement by circular permutation, i.e. sequence S′=(S₂ S₁) (FIG. 3c), correlation of S_(cyc) by S′ also gives peaks spaced out by a time T_(s) (FIG. 3d). These peaks exhibit a time lag 2^(n−1)T_(c) corresponding to the length of S₁ in relation to the peaks of the S_(cyc)*S correlation.

[0055] Land seismic prospecting—vibration time and listening time

[0056] In land seismic prospecting, the source vibrates at a fixed point and its vibration is interrupted long enough to allow displacement thereof to a next position and coupling with the ground. The relations between the length of the elementary sequence cycle, the vibration time and the listening time therefore have to be determined. What is referred to as <<listening time>> hereafter is the maximum travel time, between a source and a wave receiver, of an elastic wave that will reflect on the deepest reflector in the zone explored.

[0057] We consider an underground reflector at the time t (time required for the sound wave to move from the source to the receiver(s) by reflecting on the reflector). If the source emits a cyclic vibration S_(cyc) of elementary sequence S, this vibration will be recorded from the time t. Correlation by the elementary sequence S will give a peak at the time t, i.e. a signal between t-T_(c) and t+T_(c), where T_(c) is the period of the carrier, other peaks spaced out by multiples of the length of elementary sequence T_(s), and a minimum noise between the peaks, therefore between t−T_(s)+T_(c) and t+T_(s)−T_(c).

[0058] For a reflector at the time 0, we will thus have a narrow signal from 0 to T_(c) and a minimum noise plateau up to T_(s)−T≧T_(e). For a reflector at the time T_(c)<t≦T_(e)≦T_(s)−T_(c), we will have a narrow signal between t−T_(c) and t+T_(c) and a minimum correlation noise plateau from 0 to t−T_(c) and from t+T_(c) to T_(e). Correlation of S_(cyc) by its elementary sequence S allows to find again the series of reflection coefficients convolved by a three-arch signal of total width 2T_(c).

[0059] For a correlation over a time equal to listening time T_(e), the signal emitted can be limited to a portion of cyclic signal S_(cyc) of length T_(s)+2T_(e) equal to the sum of the length of the sequence and of twice the listening time, the elementary sequence occupying the center of this cyclic signal portion.

[0060] Similarly, it is possible to emit a signal consisting of a single elementary sequence and to correlate the record of the signals reflected by a cyclic repetition portion of this elementary sequence of length T_(s)+2T_(e), this elementary sequence occupying the center of this cyclic signal portion.

[0061] We consider the case where two sources are simultaneously recorded, one vibrating with signal S of length at least equal to twice the listening time, split up into two subsequences S1 and S2, each subsequence being equal to or greater than the listening time, the other vibrating with the sequence modified by a time lag and completed by circular permutation, consisting of the two sequences S2 and S1. The mechanism of separation of these two records is illustrated by FIGS. 4a to 4 e.

[0062] Consider a reflection coming from one of the two simultaneous sources which emits, for example, a signal S whose length is greater than twice the listening time, that can be split up into two subsequences, each one longer than the listening time, S1 and S2 (FIG. 4a). This reflection will be a signal similar to the signal emitted, weighted by the reflection coefficient and starting at arrival time t. The corresponding correlation signal (FIG. 4b) will be constructed by adding, after emission sequence S1S2, a portion of subsequence S1 starting at the beginning of the subsequence and whose length is at least equal to the listening time, and, before this central sequence, a portion of subsequence S2 whose length is at least equal to the listening time, and ending at the end of the subsequence. Correlation of the reflection by this correlation signal (S2(Te) S1 S2 S1(Te)) over a length equal to the listening time will lead to a peak at the time t (FIG. 4c) when central sequence S1S2 of the correlation signal is opposite the reflection. Correlation of the reflection coming from the first source by the correlation signal corresponding to signal S2S1 modified by a time lag, emitted by the second source, constructed by modifying the correlation signal of the first source by the length of S1 and by completing by circular permutation (FIG. 4d), will give, during the listening time, only a minimum correlation noise. Similarly, correlation of a reflection coming from the second source by the correlation sequence corresponding to the second source will lead to a correlation peak at the arrival time of this reflection, whereas correlation by the correlation sequence corresponding to the first source will only give a correlation noise maintained at a minimum level.

[0063] Correlation, over the length of the listening time, of the records where are superposed the reflections of the waves emitted by two simultaneous sources, alternately by the correlation sequences corresponding to each of the two sources thus allows to separate the reflections coming from both sources.

[0064] Generalization to m sources

[0065] The previous result can be generalized for m vibrators. Consider m vibrators vibrating simultaneously for a time T_(v)=(m−1)T_(e). The signal emitted by the first vibrator is an elementary sequence whose length is at least equal to the product of the number m of vibrators vibrating simultaneously by the listening time, i.e. T_(s)≧mT_(e). This sequence can be split up into m parts of length S₁, S₂, . . . , S_(n) greater than or equal to the listening time. Correlation of the record of a reflected signal, starting to arrive at the time t, by correlation signal (S_(n), S₁, S₂, . . . , S_(n), S₁) gives a peak at the time t and a very weak correlation noise everywhere else between 0 and T_(e). Correlation of the record between 0 and T_(e) by any one of the sequences deduced from the first sequence by circular permutation with a time lag such as a multiple of T_(e): (S₂, S₃, . . . , S_(n), S₁), (S₃, . . . , S_(n), S₁, S₂), . . . , (S_(n), S₁, S₂ . . . , S_(n-1)) will only give a minimum correlation noise.

[0066] If a second vibrator vibrates at the same time as the first one with a sequence modified by a time lag and completed by circular permutation S′=S₂ S₃ . . . S_(n), correlation by S will give only a minimum correlation noise. Only correlation by S′ will lead to a peak at the time t.

[0067] The record obtained by simultaneously recording several vibrators can thus be split up into components connected to each vibrator by successive correlations with sequences modified by a time lag equivalent to the listening time and completed on the right and on the left by the closest subsequences from the point of view of the circular permutation or, if these subsequences are greater than the listening time, by portions of a length equal to the listening time.

EXAMPLE

[0068] We consider five vibrators vibrating simultaneously with BSML sequences of 2¹⁰−1=1023 terms modulating a 60 Hz carrier and deduced from one another by time lag and circular permutation. The length of the sequence is thus $\frac{1023}{60} \sim 17$

[0069] seconds. The listening time will correspond (to within some carrier periods) to one fifth of the length of the sequence, i.e. $\frac{17}{5} \sim 3.4$

[0070] seconds. The length of the (simultaneous) vibration of the vibrators will be 17 seconds. The correlation signal will have the length as follows:

17+2×3.4=23.8 s.

[0071] Minimum correlation sequences

[0072] If we accept a correlation noise of −30 dB, it is possible to use for example an emission device comprising 4 lines of 5 vibrators each, vibrating simultaneously and controlled by a 34-s sequence, for a listening time of 6.8 s. On a first line, the vibrators vibrate simultaneously or not, with the same sequence modified by a suitable time lag as described. The next vibrator lines are controlled by sequences belonging to the set of minimum correlation sequences to which the first sequence belongs. The vibrators of each line start as soon as they are in the field and ready to vibrate. It is not necessary to synchronize the start of the various lines. An appreciable productivity gain can thus be obtained.

[0073] The implementation device (FIG. 5) comprises a control unit 1 intended to apply to m vibrators coupled with the formation phase modulated periodic signals, a system 2 intended for acquisition and recording of the signals picked up by seismic receivers R1-Rk coupled with the formation in response to the periodic signals emitted, and a processing system 3 for processing the seismic signals picked up by receivers R1 to Rk, such as a programmed computer, by correlation with the signals emitted.

[0074] Control unit 1 comprises a periodic signal generator 4 and a modulation set MA intended to generate, from the periodic signals, m phase modulated periodic signals for application to the m vibrators V1-Vm.

[0075] According to a first embodiment (FIG. 6), modulation set MA comprises means 5 for generating at least one pseudo-random control sequence, elements 6 for producing a time lag so as to form, from said sequence, at least one set of m pseudo-random binary sequences with a time lag between each one respectively. The time lags are distributed over the length of the control sequence according to the number m of vibrators of each group.

[0076] A phase modulation device 7 is connected to signal generator 4 and to elements 6 intended to apply a time lag to the periodic signals, and it generates m phase modulated periodic signals. Device 7 is connected by connection means (L) to the various vibrators (V1-Vm) of each group.

[0077] According to a second embodiment (FIG. 7), modulation set MA comprises means 5 for generating at least one pseudo-random control sequence, a modulation device 7′ for modulating the periodic signals generated by generator 4 by said sequence, means 6′ for producing m modulated periodic signals comprising a time lag between each one. Device 7′ is similarly connected to the various vibrators (V1-Vm) of each group by connection means L.

[0078] Any type of pseudo-random sequence can be selected for phase modulation of the signals. 

1) A method intended for prospecting or monitoring of an underground formation by elastic waves, comprising simultaneous emission, in the ground, of seismic waves by several seismic sources (V1 to Vm) controlled each by a signal consisting either of an elementary sequence formed by coding of a signal by pseudo-random sequences, the length of each elementary sequence being at least equal to the product of the number of seismic sources vibrating simultaneously by the listening time, or respectively of this elementary sequence extended, before and after, with parts whose length is at least equal to the listening time, reception and recording of the signals reflected by the subsoil discontinuities in response to the signals emitted, and processing of the signals recorded, the respective contribution of each one of the various seismic sources being separated by correlating the signals received and recorded either by signals constructed from said elementary sequence, extended, before and after, with parts whose length is at least equal to the listening time, obtained by circular permutation, or respectively of said elementary sequence, characterized in that the signals applied to the various vibrators are formed from a periodic signal whose amplitude and derivatives with respect to the amplitude time are equated to zero at the beginning and at the end of each period, and such that its integral over a period is zero. 2) A method as claimed in claim 1, characterized in that a signal of frequency f, of the form sin 2π.ƒ.t.(1−cos 2π.ƒ.t), is used as the periodic signal. 3) A method as claimed in claim 1, characterized in that the elementary sequences are formed by phase modulation of the periodic signal and these sequences are applied to vibrators. 4) A method as claimed in claim 3, characterized in that the same control sequence is applied to all the vibrators, with a time lag and a circular permutation, the lags between two vibrators being at least equal to the listening time, and the length of the elementary control sequence being at least equal to the sum of the lags and of the listening time. 5) A method as claimed in claim 3, characterized in that the vibrators are controlled by a connected set of several minimum crosscorrelation sequences, each vibrator emitting a vibration with its own sequence independently of the other vibrators. 6) A method as claimed in claims 4 and 5, characterized in that the vibrators are divided into several groups, and a pseudo-random sequence belonging to the same group of minimum crosscorrelation sequences is assigned to each group, the various vibrators of a group being controlled by the same elementary sequence with said time lags. 7) A device intended for prospecting or monitoring of an underground formation by elastic waves, comprising at least a group of m vibrators (V1-Vm), a control unit (1) for controlling the various vibrators by phase modulated periodic signals, seismic receivers (R1-Rk) coupled with the formation, a system (2) intended for acquisition and recording of the seismic signals reflected by the subsoil discontinuities in response to the periodic signals emitted and a system (3) for processing the seismic signals recorded, by correlation with part of the periodic signals emitted, characterized in that control unit (1) includes a signal generator (4) and a modulation set (MA) for forming m elementary sequences phase modulated by a pseudo-random sequence and connection means (L) for applying simultaneously to the various vibrators (V1-Vm) of each group the signals produced by the modulation set. 8) A device as claimed in claim 7, characterized in that modulation set (MA) comprises means (5) for generating at least one pseudo-random control sequence, time lag means (6) for forming at least one set of m elementary sequences with a time lag between each one respectively, means (7) intended for phase modulation of the periodic signals generated by said generator (4) respectively by the m elementary sequences, generating m phase modulated periodic signals. 9) A device as claimed in claim 7, characterized in that modulation set (MA) comprises means (5) for generating at least one pseudo-random control sequence, modulation means (7′) for modulating, by said sequence, the periodic signals generated by generator (4), time lag means (6′) for producing m modulated periodic signals with a time lag between each one. 10) A device as claimed in any one of claims 7 to 8, characterized in that control unit (1) is suited to emit during a time depending on the length of each control sequence and on a received vibrational signal listening time. 11) A device as claimed in any one of claims 7 to 10, characterized in that it comprises p groups of vibrators, control unit (1) being suited to generate p minimum crosscorrelation elementary sequences and to apply these sequences to the vibrators of each group with a time lag. 12) A device as claimed in any one of claims 7 to 11, characterized in that the vibrators are permanently installed for seismic monitoring of an underground zone. 